Variational cluster approach to spontaneous symmetry breaking: 
The itinerant antiferromagnet in two dimensions 
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Based on the self-energy-functional approach proposed recently [M. Potthoff, Eur. Phys. J. 
B 32, 429 (2003)], we present an extension of the cluster-perturbation theory to systems with 
spontaneously broken symmetry. Our method applies to models with local interactions and accounts 
for both short-range correlations and long-range order. Short-range correlations are accurately taken 
into account via exact diagonalization of finite clusters. Long-range order is described by variational 
optimization of a ficticious symmetry-breaking field. In comparison with related cluster methods, 
our approach is more flexible and, for a given cluster size, less demanding numerically, especially 
at zero temperature. An application of the method to the antiferromagnetic phase of the Hubbard 
model at half-filling shows good agreement with results from quantum Monte-Carlo calculations. We 
demonstrate that the variational extension of the cluster-perturbation theory is crucial to reproduce 
salient features of the single-particle spectrum. 

PACS numbers: 71.10.-w, 75.10.-h, 74.20.-z 



I. INTRODUCTION 

Several materials studied in condensed-matter physics 
display unusual physical properties which are ascribed to 
strong electron correlations. -'^'^■^ In particular, these may 
give rise to rich phase diagrams with different compet- 
ing types of short-range correlations and with different 
symmetry-broken phases. Realistic many-body models 
of these materials contain an interaction-energy term of 
the same order of magnitude as the kinetic energy or even 
larger. This implies that it is quite generally inappropri- 
ate to treat these models by conventional weak-coupling 
perturbation theory or by static mean-field decouplings 
such as the Hartree-Fock approximation. 

A complementary approach consists in an exact treat- 
ment of the interacting part while the kinetic energy is 
accounted for perturbatively. For Hubbard-type models, 
this corresponds to an expansion in powers of the hop- 
ping t around the atomic limit, and leads - at the lowest 
order - to the so-called Hubbard-I approximation.^ An 
expansion in t can be organized in a systematic diagram- 
matic seriesi^iSil This approach, however, not only fails 
for t of the order of the Hubbard repulsion J7, but also for 
low temperatures, due to the degeneracy of the ground 
state. 

An interesting extension of this strong-coupling ex- 
pansion consists in dividing the lattice into clusters of 
sufficiently small size such that they can be treated ex- 
actly, followed by an expansion in powers of the hop- 
ping between the clusters. ^•^•"'^^ The expansion in the 
inter-cluster hopping can be formally carried out up to 
arbitrary order following the diagrammatic method of 
Refs. Wm . However, going beyond the lowest order, 
is quite demanding numerically and has not been car- 
ried out so far for two-dimensional systems. For one- 
dimensional (infinitely long) chains, on the other hand, 



such an expansion could be accomplished in fact ~ to 
infinite order in powers of the inter-chain hopping. 

The lowest order of the strong-coupling expansion 
in the inter-cluster hoppingS has been termed "cluster- 
perturbation theory" (CPT)..^ Actually, the CPT can be 
considered as a systematic approach with respect to the 
cluster size, i.e. it becomes exact in the limit iVc — > oo 
where A^c is the number of sites within a cluster. From 
this point of view the CPT represents an attractive 
method which is simple conceptually but nevertheless in- 
cludes short-range correlations on the scale of the cluster 
size. Moreover, the CPT provides results for an infinitely 
extended system. Consequently, the CPT Green's func- 
tion is defined for any wave vector k in the Brillouin 
zone, contrary to common "direct" cluster calculations 
for which only a few k points are available. The nec- 
essary numerical effort is moderate: Once the Green's 
function of a cluster of a given size has been calculated 
via a numerical method, e.g. the Lanczos technique, the 
determination of the lattice Green's function is numeri- 
cally much less demanding as this requires the inversion 
of a certain number of matrices with a dimension given 
by Nc only. 

CPT results for static quantities as well as for the 
single-particle spectral function have been shown to 
agree well with different exact analytical and numeri- 
cal results. ^>i2' Recently, a generalization of the method 
with different cluster shapes has successfully been used 
for an analysis of the stripe phase in high-temperature 
superconductors.^^ On the other hand, there is also a se- 
rious disadvantage of the CPT at this level: Namely, the 
method does not contain any self-consistent procedure 
which implies that symmetry-broken phases cannot be 
studied (the case of a degenerate ground state, as e.g. for 
the Hubbard-I approximation, represents an exception). 

This generates the motivation for the present pa- 
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per. We will present an extension of the CPT which 
is based on the self-energy-functional approach (SFA) 
proposed recentlypi^ The SFA provides a general vari- 
ational scheme to use dynamical information from an ex- 
actly solvable reference system (an isolated cluster) to 
approximate the physics of a system in the thermody- 
namic limit. Using the SFA it is possible to construct a 
self-consistent or variational cluster-perturbation theory 
("V-CPT") which allows to study phases with sponta- 
neously broken symmetry. The V-CPT applies to arbi- 
trary Hubbard-type lattice models with the restriction 
that the interaction be local. 

Self-consistent cluster methods can also be con- 
structed as generalizations of the dynamical mean- 
field theory (DMFT)i^*i^ as has been shown in recent 
ygj^j.gj£Jii8j^,20^ Similar as the CPT, the cellular dy- 
namical mean- field theory (C-DMFT)-''^ is based on a 
real-space formulation. The C-DMFT performs a self- 
consistent mapping of the lattice problem onto an ef- 
fective cluster model with A'c > 1 correlated sites and 
reduces to the standard DMFT for A'^c = 1- CPT and 
C-DMFT differ with respect to the concept of "bath" 
sites. The effective cluster model which is considered in 
the C-DMFT contains an infinite number of additional 
uncorrelated ("bath") sites attached to each of the A^c 
original correlated sites in the cluster. The bath param- 
eters are determined from a self-consistency condition. 
This construction ensures an optimal description of the 
local (temporal) degrees of freedom but complicates the 
method considerably. A numerically exact evaluation of 
cluster generalizations of the DMFT could so far only be 
achieved by using quantum Monte-Carlo techniquesi^SiS^ 

Recently, it has been pointed out?- within the context 
of the self-energy-functional approach (SFA) that both a 
(variational) CPT and the C-DMFT can be considered 
as extreme limits (ng = 1 and Ug = oo) of a more general 
cluster method where reference is made to an effective 
cluster model with A"c correlated sites and Ug — 1 addi- 
tional bath sites per correlated site. Hence, the SFA for- 
mally unifies the different cluster approaches and thereby 
places our proposed method in a more general context. 
SFA-cluster calculations for the one-dimensional Hub- 
bard model24 strongly suggest that it is more efficient 
to use a cluster as large as possible and set ng = 1 (no 
bath sites) - as compared to a smaller cluster and > 1 . 
This is contrary to the opposite limit of infinite dimen- 
sions: For D — oo the exact theory (namely DMFT) is 
obtained for A^c = 1 and rig = oo. 

It is, therefore, particularly interesting to apply the V- 
CPT to the two-dimensional case and to compare with 
available numerically exact results. The low-temperature 
antiferromagnetic phase of the D = 2 Hubbard model 
at half-filling represents an optimal playground to study 
the strengths and limitations of the method. The rea- 
son is that both the effects of short-range correlations 
and long-range antiferromagnetic order manifest them- 
selves in static thermodynamic quantities as well as in 
the single-particle excitation spectrum. 



The paper is organized as follows: The variational gen- 
eralization of the CPT is introduced in the following sec- 
tion. In Sec, mil we present our results for the antiferro- 
magnetic D = 2 Hubbard model. One- and quasi-one- 
dimensional systems are discussed briefly, in addition. 
The performance of the method is analyzed by compar- 
ing with numerical results from different methods. Em- 
phasis is given to the single-particle excitation spectrum. 
Finally, our conclusions and an outlook are presented in 
Sec.HVl 



II. VARIATIONAL CPT 

Consider a system of interacting electrons on a lattice 
with a Hamiltonian H consisting of a single-particle (non- 
interacting) term Hq and an interaction term Hi. We 
require that the interacting part be local. This allows for 
a partitioning of the lattice into non-overlapping clusters 
of finite size which are not connected by Hi. In the 
simplest case. Hi describes an on-site Hubbard repulsion. 
After having divided the lattice into clusters (labeled by 
R), the Hamiltonian can be written in the form 



R 



R,R' 



where 



Ht''''\R)^Y.ta.>,C^naCR>, 



(1) 

(2) 



is the non-interacting, intra-cluster part of the Hamil- 
tonian, and Hi(R) is the intra-cluster interaction part 
(which we do not need to specify). The remaining term 



Ht"''\R,R')=J2^Ra,R'tC^naCR'b, 



(3) 



is a non-interacting part connecting different clusters 
(inter-cluster hopping). The labels a, b indicate positions 
within a cluster as well as other (spin and orbital) degrees 
of freedom. CRa annihilates an electron with quantum 
numbers a within the cluster R. For simplicity, transla- 
tional invariance with respect to the "superlattice" vector 
R is assumed. 

We are interested in the single-particle Green's func- 
tion GRa.R'b{i^) = {{cRa;c'R,b))Lo- Exploiting transla- 
tional invariance and performing a Fourier transforma- 
tion to the reciprocal space, the Green's function be- 
comes diagonal with respect to the wave vector Q from 
the (reduced) Brillouin zone corresponding to the super- 
lattice. In reciprocal space, the Green's function is a 
matrix Gq{uj) with elements GQ^a.bi^) labeled by the 
cluster variables a and b. 

Let us define a "reference system" with Hamiltonian 



H' where the inter-cluster hopping 
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is switched 
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H' describes a system of decoupled clusters of finite size. 
For not too large clusters, this system can be solved ex- 
actly, and its Green's function G'{lo) can be computed 
by conventional methods such as exact diagonalization 
(ED) or quantum Monte-Carlo (QMC). Generally, the 
corresponding Green's function G'{lj) is a matrix with 

this matrix is 



(inter) 



indices {Ra) and {R'b). As TJq 
diagonal (and constant) with respect to R: 



(5) 



Within the CPT approximation, the Green's function 
G{u!) of the full problem H is expressed in terms of G'{llj) 
and the inter-cluster hopping VRa.R'b by an RPA-type 
expressiontSiS 

GQiu;) = [G'iu;)-'-VQ]-' . 

Here, Gq{uj), G'{uj), and Vq are matrices in the cluster 
indices a and h. The Fourier-transformed inter-cluster 
hopping is given by: 



1 T/ iQ-{R-R') 

-j- VRaM'b e 

R,R' 



where L is the number of superlattice sites. 

The above formalism constitutes the "usual" CPT ap- 
proach. We like to stress that the method is based on the 
exact solution of finite-size clusters in which spontaneous 
symmetry breaking cannot occur. Furthermore, it does 
not include any self-consistent procedure. Consequently, 
symmetry-broken phases cannot be studied within the 
usual CPT. 

Our proposal for a proper generalization of the CPT 
is the following: First, one should note that in the CPT 
the perturbation term is quite arbitrary and can be taken 
as any one-particle operator. The partition of the non- 
interacting part of the Hamiltonian Eq. gives us a 
certain amount of freedom that we can exploit to seek 
for an optimized starting point. As a matter of fact. 



(intra) 




any local single- 

^(inter) 



one has the freedom to add to H^ 

particle term which is then subtracted in H^^'^"^'. In 
other words, the Hamiltonian Q is obviously invariant 
under the transformation 



(intra) 



(R) 



H^"''"'>{R)+A{R) 

^ (inter 



where A(i2) is an arbitrary intra-cluster single-particle 
operator which can be expressed as 



a,b 



Formally, the rest of the procedure remains unchanged. 
If the perturbative approach was exact, the results 

(4) would not depend on A at all. As a matter of fact, this 
can easily be seen in the non-interacting limit in which 
the CPT becomes exact4ii Here the perturbation cancels 
out. 

In the interacting case, the result does depend on 
A. However, this is not a shortcoming. On the con- 
trary, this allows us to "optimize" the results of the CPT 
calculation.'*^ Indeed, we may think of choosing A such 
that the single-particle dynamics of the cluster problem 
is "as close as possible" to the exact dynamics of the 
lattice. In this way, one can hope that the perturbative 
correction is small and that the result of the perturba- 
tive calculation is accurate enough. The question is how 
to perform this optimization in practice. Note that the 
answer to this question also solves our original problem 
as one may choose the perturbation A to represent a 
"ficticious" symmetry-breaking field term (a staggered 
magnetic field, for example) since this has the form of a 

(5) one-particle operator. 
A straightforward idea to optimize A (or the strength 

h of the symmetry-breaking field) would be to express 
a thermodynamical potential (the grand potential Q, for 
example) in terms of the CPT Green's function G, which 
depends on A in turn, and to minimize the function 
(7) ri[G(A)] with respect to A. However, the following se- 
rious problems arise: As the CPT Green's function is 
approximate and as there are different ways to obtain 
the grand potential, the procedure is not unique. So 
there are several ways the potential can depend on the 
perturbation, and the results will depend on the respec- 
tive choice. Moreover, once the grand potential is given 
in terms of A, there is no physical reason to minimize 
the grand potential as this would require a correspond- 
ing variational principle to be valid which is generally 
not the case. Below, however, we will show that an ap- 
propriate variational principle can be found in fact, and 
a corresponding potential, Eq. Hll|) . can be constructed 
the stationary point of which gives an optimized A. 

Exact variational principles of the form (5ri[G] = or 
(551 [S] = where S is the self-energy are actually known 
for a long time from standard diagrammatic theoryjS^ 
The problem is that the functional dependence 0[G] or 
fi[S] is not given explicitly but has to be constructed via 
an infinite sum of renormalized skeleton diagrams. This 
has impeded the use of the variational principles in their 
original form. 

Here, the help comes from the self-energy-functional 
approach (SFA)*^ proposed recently. The SFA provides 
a way to exactly evaluate the functional - even if the 
functional dependence is not explicit. This is achieved at 
the cost of a restriction of the domain of the functional, 
i.e. the functional ri[S] can be evaluated exactly on a cer- 
tain subspace S of trial self-energies. The idea is then to 
perform a search for the stationary point on the restricted 
(9) space S. The subspace S consists of all S which are exact 
self-energies of a reference system. Clearly, the Hamilto- 
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nian of the reference system H' must be exactly solvable 
so that one is able to compute the self-energy in practice. 
Furthermore, general argumentsi^ require that H' must 
have the same interaction part as H. The one-particle 
part of the reference system, however, is completely ar- 
bitrary and its parameters may be used to optimize the 
trial self-energy. Note that these conditions are fulfilled 
for the case considered here. Constructing the reference 
system by dividing the lattice into small clusters, the SFA 
concept just yields the desired variational CPT. 

To verify this, it is sufficient to realize the following: 
Equation ® which approximates the Green's function 
Gq{ijj) of H in terms of the Green's function G'(w) of 
the system of decoupled clusters H' and the inter-cluster 
hopping Vq, can be cast into the form of a Dyson equa- 
tion, 

GQH = (Gg)H-i-SH)-i. (10) 

Here g'^q {uj) = (w -f — t — Vq)~^ is the free Green's 
function of the infinite lattice given in terms of the chemi- 
cal potential /j,, and the intra- and inter-cluster hopping t 
and Vq, respectively, and S(w) is the cluster self-energy. 
One can therefore state that the CPT consists in approx- 
imating the self-energy of the lattice problem H by the 
(Q independent) self-energy S(ci;) of the reference system 
H'. 

The optimization problem mentioned above is now 
solved in the following way: Varying A corresponds to 
varying the one-particle parameters of the reference sys- 
tem H' , the interaction part being kept fixed. For any A, 
the reference system can be solved to get the self-energy 
S. Thus, the self-energy is parameterized as S = S(A). 
Furthermore, the Green's function G' = G'(A) and the 
grand potential f2' = ri'(A) of the reference system H' 
can be calculated. Following Ref. 0, the self-energy- 
functional for the trial self-energy S(A) can be evalu- 
ated exactly, i.e. r2[S(A)] can be calculated. This yields 
a function il(A) = r2[S(A)] the explicit form of which 
is taken from Ref. El 

r2(A) = f]'(A) 

+ ^yytrln— — 

GgHzc.„)-i-S(A,zc.„) 

- Lr^trln(-G'(A,zc^„))- (H) 

n 

Here, the frequency sums run over the discrete Mat- 
subara frequencies iujn, L is the number of clusters (or, 
equivalently, the number of Q points), and bold symbols 
denote matrices with respect to cluster indices a and h. 
Note that for the evaluation of the grand potential Hll() 
one needs the CPT Green's function, Eq. Searching 
for the stationary point of the function ri(A) means to 
search for the stationary point of the exact self-energy 
functional on the restricted domain of TJ'-representable 
self-energies. This prescription tells us which approxi- 
mate cluster self-energy as best as possible describes the 
exact one. 



III. RESULTS 

We have applied the variational CPT (V-CPT) pre- 
sented above to the single-band Hubbard model at half- 
filling and zero temperature. The Hamiltonian reads: 

i/ = ^tr,r' cj. o-c^'^c + nr,|nr,j, • (12) 

r.r' V 

Here Cr.a annihilates an electron with spin projection 
a =tji at the lattice site r, rir^a — cl ^Cr^a, and the 
hopping tr^r' is equal to t and nonvanishing for nearest- 
neighbor sites (r, r') only. Throughout the paper, t = 1 
sets the energy scale. 

In principle, the "best" result is obtained by using a 
completely general single-particle term A. However, this 
would imply the computation of for a too large 

number of parameters making the problem numerically 
impractical. For this reason, it is more convenient to start 
with a "guess" of the appropriate physical symmetry- 
breaking field. For half-filling, a good candidate is cer- 
tainly a staggered field producing a Neel ordered state. 
For simplicity, we consider clusters containing an integer 
number of antiferromagnetic unit cells. With the nota- 
tion a = (r, (t), the corresponding A has the form 

Aa,6 = h 6a.b Za Vr , (13) 

where rjr = +1 (= —1) on sites of sublattice A (B), 
Za = ±1 for spin projection tr =t, i, and h is the strength 
of the ficticious staggered field. The optimal value of h 
will be obtained by minimization"*^ of ^{K) — ri[S(/i)] as 
given by Eq. ifTT)) . Obviously, h — Q corresponds to the 
usual CPT approximation. We stress again that via the 
transformation ||HJ), the staggered field is strictly equal to 
zero in the original Hamiltonian It only appears in 
an intermediate step in the Hamiltonian of the reference 
system H' to parameterize the trial self-energy. Thus, 
/i is a variational parameter without a direct physical 
meaning in the original lattice Hamiltonian H . However, 
it does introduce a true staggered field in the reference 
(cluster) Hamiltonian H' . 

For the numerical calculations we first consider a de- 
composition of the lattice into "-s/TO x -^10" clusters as 
indicated in Fig. ^ Following Ref. 0i open boundary 
conditions are used. To evaluate the self-energy func- 
tional, the grand potential ^'[K) and the Green's func- 
tion G'{h) for a cluster are computed using the stan- 
dard Lanczos algorithm. ■^^ The self-energy is obtained as 
T,{h) = G'oih)-^ - G'{h)-\ The sum over Matsubara 
frequencies in Eq. (|ll|l can be transformed into an inte- 
gral over real frequencies.^^ After frequency integration 
and Q summation, we obtain Q{h) from Eq. Hll|l . A 
Lorentzian broadening lu uj + i6 with finite (5 = 0.1 is 
used. For this choice typically 500 Q points are sufficient 
for convergence of the results. We have checked that the 
results do not significantly depend on 6. 
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FIG. 1: Decomposition of the D — 2 square lattice into 
VTO X vTo clusters. 



A. Static quantities 

Using Eq. the grand potential Q{h) = Q[T,{h)] 

has been calculated for ?7 = 8 at half filling and T = 0. 
The result is shown in Fig.|21 As anticipated, O depends 
on h. Three stationary points are found: a maximum at 
h — and two (equivalent) minima for nonvanishing val- 
ues h « ±0.18. This means that the interacting system 
"prefers" a symmetry-broken state with a non-vanishing 
staggered magnetization m, as one would have expected 
physically. 

The stationary point at ft, = corresponds to the 
usual CPT. For h — the ground state of a single 
cluster shows antiferromagnetic correlations but is non- 
degenerate. Hence, the cluster Green's function and the 
self-energy are spin independent. This implies that there 
is no coherent continuation of the antiferromagnetic cor- 
relations across the cluster boundaries within the usual 
CPT. Consequently, the order parameter m = 0. The 




FIG. 2: Dependence of the grand potential Q (per site) on the 
ficticious staggered field h as obtained by evaluating the self- 
energy functional Q — Q['S(h)]. The lattice is decomposed 
into "vTO X vTO" clusters (see Fig. . Parameters: on-site 
repulsion U = 8, temperature T — 0, half-filling. The optimal 
staggered field is found to be /i = ±0.18. The energy unit is 
given by the nearest-neighbor hopping. 
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FIG. 3: U dependence of the ground-state energy per site 
Eo at half-filling and zero temperature as obtained by dif- 
ferent methods: direct exact diagonalization (squares), usual 
(non- variational, h — 0) CPT (diamonds), and variational {h 
optimized) CPT (V-CPT, circles) for ^TOx \/T0 clusters. Ad- 
ditionally, results from a variational Monte Carlo calculation 
(crosses)^ and a QMC simulation (error bars)— are shown 
for comparison. VMC and QMC results are extrapolated to 
Nc = (X. VMC error bars are smaller than the symbol size. 



binding energy per site that is gained by a coherent 
matching of antiferromagnetic clusters for finite h, can 
be read off from Fig. [3 to be An « 0.043. This is small 
as compared to | J| — At'^/U — 0.5 as there are contribu- 
tions from bonds connecting different clusters only. 

From the value of the grand potential at the opti- 
mal field h the ground-state energy is obtained as Eq = 

± /i(iV). We have performed calculations for different 
U. Fig. shows Eq as a function of U for the respec- 
tive optimal ficticious field h (V-CPT) and for h — 
(CPT). The results are compared with those of an exact- 
diagonalization calculation for the isolated cluster with 
A^c — 10 sites (direct ED). Furthermore, the results 
of a variational Monte-Carlo (VMC) calculation22i us- 
ing a Gutzwiller-projected symmetry-broken trial wave 
function and the results of an auxiliary-field quantum 
Monte-Carlo (QMC) studjiSS are displayed for compari- 
son. VMC and QMC data for different cluster sizes A^c 
are extrapolatedSSi^ to A"c = oo (and to T = 0, in the 
latter case). 

As compared to the ground-state energy that is ob- 
tained by diagonalization of an isolated cluster ("direct 
ED"), the (usual) CPT result represents a considerable 
improvement, as can be seen in the figure. Note that 
CPT (and V-CPT) recover the exact result in the non- 
interacting limit. The gain in binding energy is due to 
the (approximate) inclusion of the inter-cluster hopping. 
A comparison of CPT with Monte-Carlo results (VMC, 
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FIG. 4: Double occupancy d = (n^.tn^,;) calculated as 
d = d^l{^, U)/dU, potential energy per site -Epot = U d, and 
kinetic energy per site iSkin = Eo — Spot as functions of U 
obtained for the same system as in Fig. |3 via V-CPT. 
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FIG. 5: Dependence of the grand potential fl (per site) on 
the ficticious staggered field h as obtained by evaluating the 
self-energy functional fl = f2[S(/i)]. Result for the D = 1 
Hubbard model at (7 = 8, T = and half-filling. Reference 
system: decoupled set of Hubbard chains with 10 sites each. 
The inset displays Q{h) on a finer scale. 



QMC), however, still shows a sizable discrepancy. On 
the other hand, our variational CPT method perfectly 
agrees within the error bars with both Monte-Carlo re- 
sults for Eq in the entire U range. This shows that a 
proper description of long-range order is essential to get 
the ground-state energy accurately. Note, however, that 
for the ground state itself and for dynamical quantities, 
the inclusion of short-range correlations is at least equally 
important (see below). 

Let us discuss a few other static quantities. Fig. ^ 
shows the double occupancy d = {nr-\nri) as a function of 
U. The double occupancy is obtained by numerical differ- 
entiation of the grand potential d = dil{ii, U)/dU (at its 
respective minimum value). It monotonously decreases 
from the non-interacting value d — {nr^){nri) = 0.25 
and correctly tends to approach the strong-coupling limit 
d — 0. Already for U of the order of the free band width, 
a strong suppression of d is found {d « 0.052 for U = 8). 
This indicates a quick crossover from a Slater-type (itin- 
erant moments) to a Heisenberg-type antiferromagnct 
(local moments) with increasing U. The potential en- 
ergy -Epot = U d and the kinetic energy Eq — -Epot with 
Eq — n + ii{N) and ^ = U/2 are shown in addition. 
Despite the fact that local-moment formation is almost 
completed for J7 = 8, there is still a considerable kinetic 
energy Ekin ~ —0.915. This has to be attributed to the 
residual kinetic exchange. 



B. One-dimensional case 

One may ask whether or not the variational proce- 
dure always yields an antiferromagnetic state, i.e. also in 
those cases in which this is not expected physically. For 



example, an antiferromagnetic state is prohibited in one 
dimension as quantum fluctuations break up any long- 
range spin order 1^ Mean- field methods, such as Hartree 
Fock, however, often yield a Necl state also in one dimen- 
sion. In a strict mean-field theory, spatial correlations 
are neglected altogether. Due to the inclusion of short- 
range correlations, the variational CPT is clearly supe- 
rior as compared to mean-field theory. For any finite Nc, 
however, longer-range spatial correlations are neglected. 
Hence, the V-CPT may be considered as a mean-field 
approach on a length scale exceeding the cluster dimen- 
sions. We therefore expect "mild" reminiscences of typi- 
cal mean-field artifacts. 

To test this, we have performed calculations for the 
one-dimensional Hubbard model. The reference system 
consists of a decoupled set of finite Hubbard chains with 
A^c sites each. Fig. |21 shows the grand potential as a 
function of the ficticious staggered field h for U = 8. As 
one can see, the minimum of fl is given by /i = 0, i.e. 
the V-CPT predicts the system to be a paramagnet, as 
expected physically. We conclude that for this case quan- 
tum fluctuations are taken into account in a sufficient way 
to prevent the system from becoming antiferromagnetic. 

The results are not so straightforward if one consid- 
ers a one-dimensional two-leg Hubbard ladder. The ref- 
erence system consists of decoupled finite ladders with 
A^tung rungs, i.e. Nc = A^rung x 2. Results for N^ung = 2, 
A^tung = 4 and Aiung — 6 are shown in Fig. Despite 
the fact that the system is one dimensional, the calcula- 
tions predict a finite value for the staggered field and for 
the staggered magnetization. Clearly, this is an artifact 
of the remaining mean-field character on a longer length 
scale. However, we can see from Fig. that the optimal 
value of h rapidly decreases when improving the approx- 
imation, i.e. with increasing size of the clusters in the 
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FIG. 6: Q{h) as in Fig. |S] but for an infinite one-dimensional 
Hubbard ladder. The reference system consists of finite (4x2 
and 6x2) ladders. Arrows indicate the optimal value of the 
ficticious staggered field. 



reference system. This is consistent with the fact that no 



finite magnetization is expected in the N„ 



oo limit. 
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It is interesting to see what happens if one uses the 
2 ladders in order to build up a true two- 
dimensional system. The results are plotted in Fig. [7| 
In this case, the optimal value of h is much larger than 
in Fig. El and the order parameter (see below) remains 
finite and depends only weakly on the cluster size. This 
signals that for the two-dimensional system the antifer- 
romagnetic state is genuine, in contrast to D = 1. 

Fig. [7| also shows the results for the two-dimensional 
lattice using different "square" clusters, 2x2, VS x VS 
and vTO x vTO (the \/8 x \/8 cluster is obtained by dis- 
carding the rightmost sites in the first and the third line 
of the vTO X \/lO cluster shown in Fig. ^ . The compar- 
ison shows that convergence with respect to the function 
il{h) is not yet achieved for the largest cluster size con- 
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FIG. 7: fl{h) as in Fig.inibut now coupling the finite ladders 
(4x2 and 6 x 2) to a two-dimensional square lattice. Results 
using the decomposition of the square lattice into 2x2, \/8 x 
%/8 and vTO x \/T0 clusters are shown for comparison. Arrows 
indicate the optimal ficticious field. 



sidered here. In the limit of very large clusters the SFA 
becomes formally exact as the trial self-energy S is de- 
fined to be the exact self-energy of H'. In this limit, we 
expect the location of the minima (±/io) of the function 
il{h) to go to zero, or the function il{h) to become fiat 
in a region around h ^ 0. The reason is that in the infi- 
nite system a finite value for the staggered magnetization 
will already be produced by an infinitesimally small field. 
Note that for the series of 2x2, 4x2, 6x2,... clusters, 
the reference system H' does not approach the original 
two-dimensional Hubbard model H. 



C. Order parameter 

While the staggered magnetization m for the one- 
dimensional ladder system rapidly decreases with cluster 
size, m remains finite and depends only weakly on the 
cluster size in case of the two-dimensional system. Dif- 
ferences in m are found to be less than 1-2% for the dif- 
ferent cluster geometries considered in Fig. [7| A relative 
difference Am/m < 0.005 is found when comparing the 
result for the 10-site and the 8-site cluster. The staggered 
magnetization is defined as m = dil/dhcxt in the limit 
^cxt where hext is the strength of an external physical 
staggered field (not to be confused with the ficticious field 
h). Adding a respective field term to the Hamiltonian H 
and performing the derivative with respect to /icxt of the 
grand potential (at the optimal ficticious field strength 
h), yields m = (l/N,) E.(-l)'''' (Kt) - where r 

runs over sites within a cluster, A^c is the number of clus- 
ter sites, and {rir.a) = {^^/'^) j'^^duo^^Gr,r,a{}jJ+i^^)- 
This is the usual expression for the staggered magne- 
tization, but averaged over the cluster. For the two- 
dimensional Hubbard model at ?7 = 8 we find m « 0.80. 

Fig.lHlshows the local average occupation {ur^]) for the 
sites r within the -s/TO x -s/IO cluster for U = %. As any 
cluster approximation (constructed in real space), the V- 
CPT necessarily breaks the translational symmetries of 
the lattice: The approximate self-energy is obtained from 
a translationally non-invariant reference system which re- 
sults from the decomposition of the original lattice into 




FIG. 8: Local average occupation {rir,]) for [7 = 8. 
1 — {n^,t) (not plotted). 
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FIG. 9: Comparison of the staggered magnetization m as 
a function of U at half-filling obtained by different methods: 
variational CPT (V-CPT), variational Monte-Carlo (VMC>22. 
and quantum Monte-Carlo (QMC)j— see text. The arrow in- 
dicates the result m = 0.62±0.04 of a Green's-function Monte- 
Carlo stud}^ for the two-dimensional Heisenberg model. 



decoupled clusters of finite size. This implies that the 
local Green's function, which is computed from the self- 
energy using the Dyson equation, and thus the local oc- 
cupations cannot be expected to be homogeneous (within 
a sublattice). It is interesting to see, however, that this is 
not a severe drawback: Fig. |S1 shows that the variations 
of the spin-dependent local occupation and the local or- 
dered moment are very moderate within a sublattice. 

A much more inhomogeneous state with strongly vary- 
ing local occupations is obtained when coupling the fic- 
ticious field h to two sites within the cluster only. This 
variant has been considered using the \/T0 x -\/T0 and 
the VS X V8 clusters. In this case, too, a finite optimal 
value for h and antiferromagnetic long-range order are 
found (not shown). The grand potential 57 at the op- 
timal field, however, is considerably larger than in the 
usual case where h is coupled to all sites within a cluster. 
This shows that despite the artificial breaking of transla- 
tional symmetry, a homogeneous state is restored as far 
as possible. 

The U dependence of the staggered magnetization is 
plotted in Fig.j^in comparison with the VMC results^ 
(Gutzwiller-projected symmetry-broken trial wave func- 
tion) and the results of auxiliary-field quantum Monte- 
Carlo (QMC).^° Within the QMC, the order parameter is 
obtained from simulations of the static spin-spin correla- 
tion function at low temperatures. It is assumed that the 
system effectively behaves as if at T = when the ther- 
mal correlation length exceeds the cluster dimensionsi^ 
VMC and QMC data are extrapolated to = ooi22i2a 

As one can see from Fig.|51 the variational CPT yields 
a staggered magnetization which strongly disagrees with 
QMC data. In the Heisenberg hmit [/ ^ oo, the V-CPT 
seems to predict the staggered magnetization to approach 
unity. However, whether it really becomes m = 1 in this 



limit is not clear from our calculations yet. On the other 
hand, physically, one would expect a reduction of the 
staggered magnetization due to transverse spin fluctua- 
tions. In the two-dimensional Heisenberg model, several 
methods starting from the simplest spin-wave theory up 
to Monte-Carlo methods, predict a reduction of the mag- 
netization to about 60% of its saturated valueS,. On 
the other hand, the V-CPT agrees well with the results 
of the variational Monte-Carlo study, where transverse 
spin fluctuations are not fully taken into account as well. 
Furthermore, there is a very good qualitative agreement 
with respect to the size of m and the trend of m(U), 
when comparing with the results of a dynamical mean- 
field calculationp2i (as the DMFT calculation has been 
performed for the D — oo hypercubic lattice, one has to 
rescale the energies by a factor four to obtain the same 
variance of the non-interacting density of states which 
may serve as the energy unit). 

This appears to be somewhat surprising since spatial 
correlations are neglected altogether in the DMFT and 
also in the VMC calculation where local Gutzwiller pro- 
jectors are used, while the V-CPT does include the cou- 
pling to short-range correlations on the scale of the clus- 
ter size. One has to bear in mind, however, that the size 
of the order parameter is strongly affected by the cou- 
pling to long-range spin excitations. Recall that in two 
dimensions and for any finite temperature the Mermin- 
Wagner theoremSiiS^ shows that antiferromagnetic long- 
range order is destroyed due to spin waves with wave vec- 
tor q — + 0. Hence, the overestimation of the staggered 
magnetization could be ascribed to the residual mean- 
field character of the V-CPT on a length scale exceed- 
ing the size of the cluster.'''^ This view is also substan- 
tiated by our results for the one-dimensional Hubbard 
ladder which have been discussed above: To achieve a 
clear suppression of long-range order within the V-CPT, 
reference systems (finite ladders) as large as 6 x 2 have 
been required (Fig. ^ . This is an indication that in two 
dimensions a \/T0 x -\/T0 cluster might be to too small to 
include non-negligible effects of spin excitations on the 
order parameter. 

There is another important point which has to be taken 
into account in this context: For a cluster of a given size, 
an optimal V-CPT calculation should not only consider 
the ficticious staggered field h but any one-particle term 
in the Hamiltonian of the reference system as a varia- 
tional parameter. It is in fact reasonable to assume that 
there is room for improvement: Consider, for example, 
the hopping between nearest neighbors within the cluster 
as an additional variational parameter. Actually, this has 
already been considered in Ref. for the D — 1 Hub- 
bard model. There it was found that the optimal intra- 
cluster hopping is increased as compared to the nearest- 
neighbor hopping in the original lattice although the ef- 
fect turned out to be rather weak. Here, the situation is 
different due to the antiferromagnetic long-range order: 
In the limit [/ — > oo, an increased intra-cluster hopping 
implies an increased effective exchange interaction |J|. 
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Assuming the optimal ficticious field h to be unchanged, 
this tends to decrease the order parameter m. We have 
performed corresponding calculations which show that a 
variational adjustment of the intra-cluster hopping at a 
considerably increased value is very likely in fact. How- 
ever, a conclusive result has not yet been obtained. The 
reason is that the grand potential as a function of two 
variational parameters, the ficticious staggered field and 
the intra-cluster hopping, becomes rather flat in a wide 
region around the stationary point. This requires an im- 
proved accuracy in the evaluation of the self-energy func- 
tional which is difficult to achieve with a finite value for 
the (Lorentzian) broadening parameter d. Work into this 
direction is in progress i2i 

D. Dynamical quantities 

While the V-CPT must be considered as mean-field- 
like on a length scale exceeding the cluster size, it does 
account for short-range spatial correlations as the cluster 
problem is solved exactly. For the two-dimensional Hub- 
bard model at half-filling, short-range spin correlations 
are known to manifest themselves in dynamical quanti- 
ties such as the local density of states. 

Fig.llOlshows the spin-dependent local density of states 
(DOS) Pai^^) for [/ = 8 which is calculated as a staggered 
average over the sites in a cluster: 

where prai^) = (— l/7r)ImG',.^r,(T('jJ + *0+). Roughly, the 
spectrum consists of two broad peaks around cu = ±5 
and two strong and narrow peaks at about lo — ±3. For 
both the high- and the low-energy excitations a strong 
spin polarization corresponding to m w 0.80 is clearly 
visible. There is also some finite but low spectral weight 
within the insulating gap which, however, is an artifact 
of the finite Lorentzian broadening {S — 0.1). 

The high-energy excitations in Fig. 1101 are interpreted 
as charge excitations (Hubbard "bands"). While these 
are due to local correlations, the low-energy features (at 
uj — ±3) result from (short-range) non-local correlations. 
The latter will be identified as being due to the coher- 
ent propagation of a quasiparticle, namely a "spin bag" . 
Physically, this spin bag originates from the frustration 
induced by the motion of the additional bare hole (elec- 
tron) in the antiferromagnetic spin background. The 
different spectral features can easily be identified: The 
high-energy features are due to the bare particle "rat- 
tling around" within the spin bag. This gives rise to an 
incoherent motion and broad energy "bands" , i.e. the in- 
coherent lower and upper Hubbard "bands" with a width 
set by the energy scale of the bare band width W — 8. 
The low-energy features, on the other hand, correspond 
to the above-mentioned coherent motion of the spin bag 
resulting in a strongly renormalized quasi-particle band 
with a width essentially given by 2\J\ — 8t'^/U = 1. 



Note that these peaks are absent in a mean-field ap- 
proach where off-site correlations are neglected alto- 
gether: A recent DMFT study'^* of antiferromagnetic 
order shows a rather featureless DOS consisting of the 
two (polarized) Hubbard bands only. Contrary, the ef- 
fect of antiferromagnetic short-range correlations can be 
included in a cluster extension of the DMFT. Additional 
structures appear in the DOS within the dynamical clus- 
ter approximation (DCA), for example. Some indications 
of the mentioned low-energy features can be found by 
using the non-crossing approximation (NCA) to evaluate 
the DCA.^^ For a conclusive interpretation, however, the 
effects are too weak - probably due to the limited cluster 
size (a 2 X 2 cluster in reciprocal space) and the finite 
temperatures considered. 

More elucidating is a comparison of the fc-resolved 
spectral density with available results from QMC simula- 
tions for isolated but larger clusters. In order to illustrate 
this point, we have plotted in Fig. 1111 the spectral func- 
tion Ak{uj) for U = 8 along high-symmetry directions in 
the Brillouin zone of the chemical lattice. The result is 
compared with the result from the usual CPT {h — 0) 
and with numerically exact QMC data from Grober et 
ali2i which are available for an TVc = 8 x 8 isolated clus- 
ter and finite but low temperature (T = 0.1). The spec- 
tral function Ak{uj) obtained from the maximum-entropy 
method (see Ref. jSTj) is shown in Fig. ^2 (bottom) . Since 
the spin-spin correlation length at T = 0.1 considerably 
exceeds the cluster dimensions, the QMC result can be 
considered as a good approximation to the T = limit. 
At half-filling the spectrum almost exactly respects the 
constraint Akito) = Ak+q{—Lo) with q = (tTjTt) which is 
predetermined by particle-hole symmetry. This must be 
considered as a strong check of the numerics. As for the 
finite system there is no spontaneous symmetry break- 
ing, the spectrum is spin independent and shows perfect 
translational symmetry with respect to the chemical lat- 
tice (periodic boundary conditions have been used). 

This must be kept in mind when comparing with the 
V-CPT. In the V-CPT the real-space spectral function 
AR,r.R'r',<T{uj) = -(l/7r)Im((cH,r,o-;cjfj, ,,,^))L'^'''' is spin 
dependent, and translational symmetry holds with re- 
spect to the superlattice vectors R only. For a proper 
comparison with the QMC data we therefore compute 

R,R' r,r' 

(15) 

see Fig. ^2 (middle) . If there was no antiferromagnetic 
order and no artificial breaking of translational symmetry 
due to the cluster approximation, i.e. if AR^r,R'r' .aii^) de- 
pended on the difference R + r — R' — r' only, this would 
correspond to the usual Fourier transformation. Here, 
Ak{Lu) is actually the diagonal element Afc,fc(w) obtained 
from two independent Fourier transformations with re- 
spect to both, R + r and R' + r'. Taking the diagonal 
element, ensures a positive definite result, Afc(ci;) > 0, 
and implies a spatial average over the cluster sites (see 
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FIG. 10: Spin-dependent local density of states (DOS) pa{ijj) 
(staggered average over the sites in a cluster) obtained for the 
same system as in Fig. Elvia V-CPT for (7 = 8. 



Ref . Il0|) . Due to this spatial average, the spectral func- 
tion is spin independent - even in the symmetry-broken 
state (an integer number of antiferromagnetic unit cells 
are included in a single cluster). 

The CPT spectral function is calculated accordingly 
but for h — Q (Fig^] top). This means that any sig- 
natures of long-range order are switched off in the spec- 
trum, and only short-range correlations (up to the clus- 
ter boundaries) are retained. Both, the CPT and the 
V-CPT result, respect the condition Ak{uj) = Ak+q{—ijj) 
with q = (tt, tt) due to particle-hole symmetry. Note that 
in both cases the spectral function is defined for any k 
point in the Brillouin zone, contrary to the "direct" clus- 
ter method (QMC). 

As already noted in the discussion of the local DOS, 
the V-CPT spectrum clearly consists of four spectral fea- 
tures, two high-energy "bands" which show strong damp- 
ing effects (incoherent Hubbard bands) and two narrow 
low-energy bands which represent the coherent dispersion 
of a quasi-particle (spin bag). This four-band structure 
is also evident in the QMC result. Comparing the en- 
ergetic positions, dispersions, weights and widths of the 
four spectral features, one can state that the agreement 
with the QMC spectrum is almost perfect. 

Roughly, the CPT and the V-CPT spectra appear to 
be similar but looking at finer structures it is obvious 
that the CPT predicts a spectral function which is quite 
different: First, and most important, there is no coherent 
low-energy band in the CPT spectrum. This shows up 
when comparing with the V-CPT around F for uj < Q 
(or around M for w > 0), for example: In agreement 
with the QMC result, the V-CPT predicts a dispersive 
low-energy band which extends continuously with spec- 
tral weight from F to X and which is clearly separated 
from the more incoherent feature at higher energies. On 
the other hand, in the CPT spectrum this is missing. In 
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FIG. 11: Density plot of the spectral function for the D = 2 
Hubbard model at (7 = 8, half-filling and T = as obtained 
by the CPT with ft = (top) and by the variational CPT 
with optimal ficticious staggered field h ^ Q (middle). The 
lattice is covered by vTO x vTO clusters. Bottom: QMC 
(maximum entropy) result, taken from Ref. 37, for the same 
parameters but for a finite low temperature T = 0.1 and an 
isolated 8x8 cluster. Dark (light) areas correspond to large 
(small) spectral weight. 



the T — M direction the low-energy features turn out to 
be too broad and are discontinuously split into several 
branches in the CPT spectrum. The dispersion around 
X is at variance with the QMC data. Finally, at higher 
excitation energies, several weak and almost dispersion- 
less bands can be found in the CPT spectrum while in 
the V-CPT there is a comparatively smooth incoherent 
background. We conclude that the variational procedure 
is crucial to achieve a qualitatively correct reproduction 
of the one-particle excitation spectrum and of the coher- 
ent quasi-particle band in particular. 

The physical reason is as follows. From previous QMC 
studies'^^ it is well known that the quasi-particle band 
is the dispersion of a spin bag, i.e. an additional hole 
(electron) which is dressed by the local distortions of the 
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spin order that are produced by the motion of the hole 
in the antiferromagnetic background. Since the hnear 
extension of the spin bag is about 3-4 sites only, this pic- 
ture is already captured by an exact diagonalization of 
an isolated small cluster. However, the emergence of a 
coherent band requires more, namely a coherent motion 
of the spin bag on a larger length scale. This is cap- 
tured in the QMC results for a large cluster of 8 x 8 sites. 
Of course, the perturbative treatment of the inter-cluster 
hopping within the CPT framework carries out a part of 
the job. This results in a string dispersion in the V-CPT 
spectrum with a bandwidth of about 2| J| = jU = 1 as 
can be read off from Fig. ^2 Also for the plain CPT the 
perturbative coupling of the clusters works into the right 
direction: Although the spectrum more or less consists 
of a two-band structure, there is a tendency towards the 
formation of a gap within each of the two bands, i.e. a 
coherent band tends to split off. Within the plain CPT, 
however, the motion of the dressed hole cannot be com- 
pletely coherent as there is no definite alignment of spins 
across the cluster boundary. Upon reaching the cluster 
edge, the spin bag encounters a misaligned spin with 50% 
probability and is partly reflected back inside the clus- 
ter. This partial loss of coherence explains the several 
bands at higher energies in the CPT spectrum which are 
absent in the V-CPT. The variational generalization of 
the CPT cures this problem by ordering spins antifer- 
romagnetically with help of the ficticious staggered field 
not only within but also across the cluster boundaries 
thereby allowing the coherent spin-bag propagation. 



IV. CONCLUSIONS AND OUTLOOK 

Correlated electron systems in two dimensions are in 
many cases characterized by strong local and non-local 
but still short-ranged correlations on the one hand and 
by long-range, e.g. magnetic, order on the other. Here, 
we have presented a variational extension of the cluster- 
perturbation theory which combines the exact diagonal- 
ization of isolated small clusters with a mean-field con- 
cept to build up an infinite lattice. Conceptually, the 
method is based on the recently proposed self-energy- 
functional approach (SFA) which sets up a very general 
variational scheme to use dynamical information from an 
exactly solvable reference system W (the isolated clus- 
ter) to approximate the physics of a system H (the D = 2 
Hubbard model) in the thermodynamic limit. 

We have apphed the variational CPT (V-CPT) to the 
Hubbard model at half-filling to study the antiferromag- 
netic phase at zero temperature. The diagonalization of 
Hubbard clusters of finite size (typically = 10) is per- 
formed using the standard Lanczos algorithm. In com- 
parison with results from variational Monte-Carlo and 
quantum Monte-Carlo studies, the V-CPT predicts the 
ground-state energy and related static quantities with 
high accuracy. While long-range antiferromagnetic or- 
der is obtained for the D = 2 model, the V-CPT yields 



a paramagnetic state for D = 1. This indicates that 
quantum spin fluctuations which inhibit an ordered phase 
in the D = I case are included properly. For one- 
dimensional Hubbard ladders the method in principle in- 
correctly predicts antiferromagnetic order; however, the 
staggered magnetization is small and tends to vanish 
when increasing the number of rungs in the cluster (up 
to 2 X 6). The finite but small m for the ladder system 
should be considered as a mild reminiscence of a typical 
mean-field artifact which shows up because longer-range 
spin correlations exceeding the cluster dimensions are ne- 
glected. A similar effect is seen for the D = 2 system: 
Here the approximation is even stronger because the lin- 
ear dimension of the cluster must be reduced even more 
to keep the number of cluster sites Nc constant. For the 
D = 2 system, antiferromagnetic order is expected phys- 
ically and is also found by the calculations. However, 
longer-ranged spin correlations give rise to a consider- 
able reduction of the order parameter which is not seen 
in the V-CPT for the maximum cluster size that has been 
considered. 

An important advantage of our method is that lo- 
cal and off-site short-range correlations within the or- 
dered phase can be treated exactly. This shows up 
when looking at dynamical quantities, such as the spin- 
dependent local density of states or the spectral function 
Ak (uj) ■ The spectral function as calculated from our self- 
consistent cluster approach agrees extremely well with 
the QMC (maximum-entropy) result for an 8 x 8 Hub- 
bard lattice at finite but low temperatures. In particular, 
it is possible to reproduce the dispersions, widths and 
weights of the different spectral features. This is due to 
the fact that the typical four-band structure arises not 
only from local correlations which arc captured in dy- 
namical mean-field theory, for example, but also from a 
strong coupling to off-site spin correlations. The forma- 
tion of a spin-bag quasi-particle as a hole which is dressed 
by the distortions of the antiferromagnetic spin structure 
that are introduced due its motion, is contained in the ex- 
act treatment of the cluster. In addition, the mean-field 
coupling of the individual clusters mediates the informa- 
tion on the spin order across the cluster boundaries and 
thereby gives a qualitatively correct description of the co- 
herent propagation of the quasiparticle. This is essential 
to reproduce the low-energy quasi-particle band in the 
spectrum. 

There are different interesting routes to be explored 
in the future: A straightforward but technically ambi- 
tious idea is to employ a quantum Monte-Carlo technique 
within the self-energy-functional framework. This would 
offer the possibility to use a decomposition of the lattice 
into larger clusters which is expected to be important 
for a reliable estimate of the order parameter for exam- 
ple. Another straightforward extension of the method 
concerns the number of variational parameters. For our 
present purposes the consideration of a single parameter, 
the ficticious staggered field h, has been sufficient. It is 
obvious, however, that the results should improve when 
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treating any one-particle parameter within the cluster as 
a variational parameter. This might also affect the mag- 
netic properties because of the direct link between hop- 
ping parameters and the effective exchange interactions. 
Here an improved numerical evaluation of the theory is 
required which is free from the use of broadening param- 
eters, see Refs. Il3l27l 

An important conceptual problem of the method con- 
sists in the fact that a suitable reference system can be 
found in case of local (on-site) interaction terms only. 
The reason for this restriction is that for models includ- 
ing e.g. a nearest-neighbor Coulomb interaction, the par- 
titioning of the infinite lattice into decoupled clusters is 
impossible without cutting the inter-cluster interaction. 
General arguments, ^'^ however, require that the interac- 
tion part of the Hamiltonian of the reference system be 
unchanged. Consequently, there is a need for an exten- 
sion of the theory to models with non-local interaction 
terms. This will be the subject of a forthcoming paper. ■^^ 

The V-CPT compares with the recent cluster exten- 
sions of the DMFT>iSJIiiS*i2iS2iSi Characteristic to both 
approaches is the combination of a numerically exact 
treatment of an isolated cluster with an approximate 
mean-field treatment of the coupling between different 
clusters. The V-CPT is, however, much easier to imple- 
ment numerically and in principle allows for a diagonal- 
ization of larger clusters which facilitates a proper finite- 
size scaling. The reason for this conceptual simplicity 
is the fact that there is no coupling to bath degrees of 
freedom which is inherent to the cluster extensions of the 
DMFT. Recent results^^ for the D = 1 Hubbard model 
have shown that it is in fact more efficient to consider 
larger clusters instead of a coupling to bath sites. For the 
two-dimensional case, however, there is no answer to this 
question at present, and it is unclear whether or not bath 
degrees of freedom efficiently speed up the convergence 
to the exact solution with increasing cluster size. Note 
that a hopping to uncorrelated bath sites is nothing but 
a modified one-particle part of the Hamiltonian of the 



reference system and thus completely in line with the 
general concept of the self-energy-functional approach. 
In particular, it can be very interesting to consider an 
extension of the V-CPT with a coupling to a few bath 
sites as has been done in previous studies. ^2*^ There are 
at least two reasons for that: First, this bridges the gap 
to the cellular dynamical mean-field theory (C-DMFT)iS 
which is obtained in the limit of infinite number of bath 
degrees of freedom. Second, bath sites are expected to 
become important for the study of doped systems as they 
can serve as a particle reservoir. 

For doped systems, the variational optimization of one- 
particle parameters should include the on-site energy of 
the cluster sites as well as the additional consideration 
of a few bath sites. This is necessary to realize fillings 
which are not commensurate with the cluster size and to 
achieve smooth doping dependencies in the entire doping 
range. 

In the present paper we have focused on a ficticious 
staggered magnetic field as a variational parameter to 
describe antiferromagnetic long-range order. Moreover, 
one can envisage to introduce in the same way a ficticious 
pairing field in order to study off-diagonal long-range or- 
der and superconductivity. While this introduces the nu- 
merical difficulty of diagonalizing clusters without fixing 
the particle number, such an approach offers the exciting 
perspective of analyzing the effects of short-range corre- 
lations in the superconducting phase. 
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